The effect of Zn-Sn pair-doping on M-type Sr-hexaferrite 
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We study the site occupancy and magnetic properties of Zn-Sn substituted M-type Sr-hexaferrite 
SrFei2- :E (Zno.5Sno.5) : EOi9 with x — 1 using first-principles total-energy calculations. We find that 
in a ground-state configuration Zn-Sn ions preferentially occupy 4/i and 4/2 sites unlike the model 
previously suggested by Ghasemi et al. [J. Appl. Phys, 107, 09A734 (2010)], where Zn-Sn ions 
occupy 2b and 4/2 sites. Our model predicts a rapid increase in saturation magnetic moment (M s ) 
as well as decrease in magnetic anisotropy compared to the pure M-type Sr-hexaferrite consistent 
with experimental observations. 



I. INTRODUCTION 

Strontium hexaferrite SrFei20i9 (SrM)is widely used 
as a material for permanent magnets along with other 
M-type hexaferrites MFei 2 0i 9 (M = Sr, Ba, Pb) due to 
their high Curie temperatures, large saturation magneti- 
zations (M s ), high value of coercivity, excellent chemical 
stability and low cost. Recently SrM is investigated for 
the use in magnetic recording mediap^ but due to its 
high intrinsic coercivity, pure SrM is difficult to be used 
as magnetic recording materials!^ Magnetism in SrM re- 
sults from Fe 3+ ions occupying five different interstitial 
sites in the ten oxygen layers of the crystallographic unit 
cell: three octahedral sites (12fc, 4/ 2 , 2a), one tetrahcdral 
site (4/i) and one trigonal bipyramid site (2b). Investi- 
gations to improve the magnetic properties of M-type 
hexaferrites have been done through substitution of vari- 
ous impurities JinlSJ Substitution of Fe 3+ ions with transi- 
tion metal ions have proven to be successful in improving 
the magnetic properties of SrM. Pair-doping of Fe atoms 
by Zn-NtP, Zr-CcP, Zn-SrP, Sn-M#l and Er-NP im- 
proved magnetic properties of SrM. Zr-Cd substituted 
SrM (SrFei2_ a: (Zro.5Cdo.5) a ;Oig) showed increase of M s 
only up to x = 0.2. Coercivity decreased continuously 
with increase of Zr-Cd concentration^!] while M s and co- 
ercivity steadily increased in Er-Ni substituted SrM with 
Er-Ni concentration. Substitution by Zn-Nb 18 , Zn-SrP^ni 
and Sn-Mg 2 J increased the saturation magnetization (M s ) 
and decreased coercivity. 

Using the density functional theory (DFT) Fang et 



alps investigated the electronic structure of SrM. Novak 
et al. calculated the electronic structure and exchange in- 
teractions of barium hexaferrite using DFT Pin spite of 
the importance of substituted SrM, only a few theoreti- 
cal investigations have been done on La substitution! 23 ! 24 ! 
To our knowledge, no electronic structure calculation has 
been done on experimental investigations of pair-doped 
SrM. The knowledge of site occupancy of the substituted 
atoms and the change in atomic magnetic moments in 
substituted SrM would be advantageous in understand- 
ing the effect of pair-doping on magnetic properties of 
SrM. 



In this work we use first-principles total-energy cal- 
culations to study the site preference and magnetic 
properties of Zn-Sn substituted M-type Sr-hexaferrite 
SrFei2- 2: (Zno.5Sno.5) 2 ;Oi9 with x = 1. Based on total en- 
ergy calculations, we conclude that in the ground-state 
configuration Zn-Sn ions preferentially occupy 4/i and 
4/2 sites in co ntras t to the model suggested by Ghasemi 
and co-worker d 7 * 20 !, where Zn-Sn ions occupy 2b and 4/ 2 
sites. We further show that our model predicts increase 
in saturation magnetization (M s ) as well as decrease in 
magnetic anisotropy energy (MAE) compared to the pure 
M-type Sr-hexaferrite (x = 0) consistent with experimen- 
tal observations. 
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II. METHODOLOGY 

To determine the site preference of Zn and Sn atoms 
in Sr-hexaferrite, we use first-principles total-energy cal- 
culations for configurations of SrFe^-a^Zno.sSno^xOig 
at x — and x = 1. x — configuration corresponds 
to the pure SrM without any substitution. We used a 
unit cell which contains two formula units of SrM in 
our calculations, x — 1 configuration was constructed 
by substituting Zn and Sn, one atom each, into the Fe 
sublattices of SrM. Total energies and forces were cal- 
culated using density-functional theory (DFT) combined 
with projector-aug mented - wave (PAW) potentials as im- 
plemented in VASP. 1251221 All calculations were spin po- 
larized according to Gorter's ferrimagnetic ordering of Fe 
spins J22l Brillouin zone sampling was performed using the 
Monkhorst-Pack scheme 2 ^ with a Fermi- level smearing of 
0.2 eV applied through the Methfessel-Paxton method 2 ^ 
for relaxations and the tetrahedron methocP^ for static 
calculations. Using the quasi-Newton algorithm,^ we 
performed geometric optimizations to relax the posi- 
tions of ions. Electron exchange and correlation were 
treated with generalized gradient approximation (GGA) 
as parameterized by the Perdew-Burke-Ernzerhof (PBE) 
scheme.^ To improve the description of Fe 3d elec- 
trons hybrid functionals in the Heyd-Scuseria-Ernzerhof 
(HSE) schemtPHaUand GGA+U method were employed. 
Formula by Vosko et alP^ was used for the interpo- 
lation of the correlation part of the electron exchange 
and correlation potential. A plane-wave energy cutoff 
of 520 eV was used, and volumetric optimizations were 
performed to obtain ground state of SrFe^Oig. For 
SrFei2-x(Zno.5Sno.5)a;Oi9, with x = 1, only atomic re- 
laxations were performed with a plane-wave energy cut- 
off of 400 eV. In each case, geometric optimizations were 
performed until the energy difference between two suc- 
cessive ionic optimization steps was less than 0.001 meV. 



A. GGA+U 

Although GGA is adequate for most calculations, when 
applied to systems such as transition metal oxides that 
contain 3d electrons, it underestimates the band gap 
and the magnetic moment. Therefore to improve the 
description of the Fe 3c? electrons, we used the simpli- 
fied rotationally invariant approach to the L(S)DA+U 
method described by Dudarev et al.p^with GGA instead 
of the L(S)DA exchange and correlation functional. The 
method requires an effective U value (U c s) equal to the 
difference between the Hubbard parameter U and the ex- 
change parameter J, which is part of the added penalty 
functional that represents the double counting correction 
to the total energy of L(S)DA (in our case GGA). The 
U and J parameters can either be determined from ex- 
periments or estimated using restricted L(S)DA (GGA) 
calculations. Novak et al. performed GGA+U calcu- 
lations using full potential linearized augmented plane 



wave (FPLAPW) methodP^Their calculations compared 
two values, U = 3.4 eV and U = 6.9 eV, with GGA 
which corresponds to U = eV. We used the magnetic 
moment of Fe atoms predicted by calculations using hy- 
brid density functionals to determine the U value. Using 
the determined U, the electronic structure of SrFe^Oig 
in Gorter's spin structure was calculated and volumetric 
optimizations were performed to obtain the respective 
ground states. 

B. Hybrid functionals 

Hybrid density functionals have much more accuracy 
than GGA functionals. But they are computationally 
expensive for large systems. Available Hybrid function- 
als include a certain amount of the Hartree-Fock ex- 
change. The HSE scheme is a screened Coulomb func- 
tional that is based on the GGA (PBE) exchange func- 
tional with a short-range modification. The accuracy 
of the HSE hybrid functional is comparable to B3LYP 
and PBE0 functionals whereas it considerably reduces 
the computational effort. It has been successfully ap- 
plied for large molecules and periodic systems.^ To de- 
termine the U value we calculated the electronic struc- 
ture of SrFei 2 Oigusing HSE functionals. HSE calcula- 
tion was started from the results of the GGA calcula- 
tion. The volumetric relaxed structure and charge densi- 
ties of SrFe^Oigfrom the GGA calculation were used to 
obtain the energy and magnetic moments through HSE 
calculations. The U value which reproduced the average 
magnetic moments of each symmetry site of the unit cell 
closest to the HSE calculation was then chosen. 



III. RESULTS 

A. Pure SrFei 2 Oi 9 

Pure SrM (SrFei2-a;(Zno.5Sno. s^Oig with x = 0) has 
the hexagonal crystal structure and the space group 
P63 / mmc as shown in Fig. [lj As shown in Table [TJ the 
unit cell has 11 inequivalent sites: one Sr site of multi- 
plicity 2, five Fe sites of multiplicity 2, 2, 4, 4 and 12, 
and five oxygen sites of multiplicity 4, 4, 6, 12, and 12. 
In total the double formula unit cell has 64 atomic sites. 
First, using GGA functional form wc optimized the crys- 
tal structure of pure SrM using volumetric and ionic opti- 
mizations. The optimized structure was then used to de- 
termine magnetic properties of pure SrM crystal through 
a static calculation using HSE hybrid functional. Finally 
the variation of magnetic properties with respect to effec- 
tive U parameter of the GGA+U method was determined 
through fully optimized calculations. The dependence of 
the total magnetc moment of the unit cell and the atomic 
magnetic moments of individual symmetry sites for Fe +3 
ions on the parameter U is plotted in Fig. [2] and Fig. [3] 
The values from HSE functional and the nominal value 
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FIG. 1. (Color online) One double formula unit cell of pure 
SrFei2 0ig in Gorter's spin configuration. Black spheres are 
Fe atoms, blue spheres are O atoms, and orange spheres are 
Sr atoms. 



o-o GGA+U 
--- HSE 
— Nominal 



TABLE I. Average magnetic moments of atoms in SrFei20ig 
calculated using GGA, GGA+U (with U = 3.7 eV), and HSE 
functionals. Total magnetic moment of the double unit cell is 
also given. All moments are given in /is. 



site 


GGA 


GGA+U 


HSE 


Sr(2d) 


-0.005 


-0.003 


-0.008 




O. 1 oo 




4.149 


Fe(26) 


3.541 


4.057 


4.044 


Fe(4/0 


-3.426 


-4.038 


-4.042 


Fe(4/ 2 ) 


-3.166 


-4.095 


-4.093 


Fe(12fc) 


3.719 


4.170 


4.138 


0(4e) 


0.383 


0.353 


0.391 


0(4/) 


0.121 


0.089 


0.093 


0(6h) 


0.085 


0.027 


0.031 


0(12fc) 


0.093 


0.084 


0.082 


0(12fc) 


0.175 


0.180 


0.196 


Total 


39.89 


40.00 


39.95 



calculations. 

Through volumetric and ionic optimizations, we ob- 
tained a — 5.93 A and c = 23.20 A as equilibrium lattice 
constants in good agreement with experimental values of 
5.88 A and 23.04 A, respectively.^ The total densities 
of states (DOS) for the optimized Gorter's configuration 
calculated with HSE, GGA+U (U = 3.7 eV) and GGA 
are shown in Fig. |4j HSE and GGA+U predict an in- 
sulating state while GGA predicts a metallic state with 
nonzero DOS at Fermi energy. The magnetic moments 
of individual atoms of the 11 inequivalent atomic sites 
and the saturation magnetic moment of the unit cell, 
from GGA, GGA+U and HSE functional forms are given 
m Table |U Individual magnetic moments of the Fe sites 
predicted by GGA are much less than those by HSE. 
GGA+U's prediction of the atomic magnetic moments 
is in good agreement with HSE results. Total magnetic 
moment values are similar to the nominal value of 40 \xb 
per unit cell for all three methods. 



U Parameter 

FIG. 2. (Color online) Variation of the total magnetic mo- 
ment of the unit cell with respect to the effective U parameter 
for the Goter's configuration of SrM. The values from calcula- 
tions using HSE functional and the nominal value are shown 
as horizontal lines. 

(40 [Ib) are also shown as horizontal lines for comparison. 
Fig. [2] shows that the total magnetic moment M s does 
not serve as a good criteria for determining the appro- 
priate U value as it is nearly constant over a wide range 
of U values. On the other hand, Fig. [3] shows that the 
magnetic moments on inequivalent sites vary smoothly 
over the range of U values we considered and, more im- 
portantly, they cross their corresponding values of HSE 
calculation simultaneously near U = 3.7 eV. Therefore, 
we set U to 3.7 eV for GGA+U method for all subsequent 



B. SrFei2- :E (Zno.5Sno.5) : i;Oi9 with x = 1 

For x — 1 concentration, one Zn and one Sn atom 
are substituted to two of the 24 Fe sites of unit cell 
as shown in Fig. [5] All Fe sublattices listed in Ta- 
ble [I] have more than one equivalent atomic site. In 
general substitution of Zn-Sn atoms breaks the symme- 
try among the equivalent sites of pure SrM. Out of all 
552 (= 24 x 23) possible Zn-Sn-substituted SrM struc- 
tures, the 24 symmetry operations of the space group 
P63 /mmc render substantial number of structures equiv- 
alent leaving only 55 inequivalent structures. We label 
these inequivalent configurations using the convention of 
Zn(site for Zn)Sn(site for Sn). (unique index). For ex- 
ample, when Zn and Sn are substituted to 2b and 12k 
sites, respectively, there are 24 (= 12 x 2) possible struc- 
tures but only two of them are inequivalent configura- 
tions due to symmetry. Therefore, they are labeled as 



4 



O-O 2a GGA+U 
O-O 2b GGA+U 
□ 12k GGA+U 

— 2a HSE 

— 2b HSE 
12k HSE 
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FIG. 3. (Color online) Variation of atomic magnetic moments with respect to the effective U parameter for the Goter's 
configuration of SrM ; (a) positive magnetic moments for 2a, 2b, and 12k sites and (b) negative magnetic moments for 4/i and 
4/2 sites. The values from calculations using HSE functional are shown as horizontal lines. 
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FIG. 4. (Color online) Total densities of states (DOS) 
(states/eV) of SrFei 2 Oi 9 calculated using GGA, GGA+U 
(U = 3.7 eV), and HSE methods. Positive values in red cor- 
respond to spin-up states and negative values in black corre- 
spond to spin-down states. Energy band gaps are measured 
to be 0.0, 0.93, and 1.19 eV, respectively. 



Zn(26)Sn(12/fc).l and Zn(26)Sn(12fc).2. These structures 
were constructed by substituting Zn and Sn atoms to the 
respective sites of a SrM unit cell fully optimized using 
GGA+U (17 = 3.7 eV). After substitution, the atomic 
positions of the structures and spin polarized electron 
charge densities were fully optimized and their total en- 
ergies were compared. All configurations maintained the 
original hexagonal structure. 

In Table [XT] we lists the ten lowest energy configu- 
rations among all incquivalent Zn-Sn-substituted SrM 
structures we considered. The ground state configura- 
tion Zn(4/i)Sn(4/2).l that has Zn and Sn atoms sub- 
stituted to 4fi and 4/ 2 sites, respectively, has an in- 
crease in M s of approximately 10 /is with respect to pure 





(a) (b) 

FIG. 5. (Color online) The ground state structure of Zn-Sn- 
substituted SrM (configuration Zn(4/i)Sn(4/2).l) with spins 
oriented in (a) easy axis (001) and (b) hard axis (100). Black 
spheres are Fe atoms, blue spheres are O atoms, orange 
spheres are Sr atoms, green sphere is Zn atom and light blue 
sphere is Sn atom. 



SrM consistent with the experimental results by Ghasemi 
and co-workers!^^ The structure suggested by Ghasemi 
with Zn 2+ in 2b site and Sn 4+ in 4/2 site (configura- 
tion Zn(26)Sn(4/2).l in Table |ll|) is an excited state con- 
figuration ranked #8 with a much higher total energy 
(0.277 eV). Furthermore, configuration Zn(26)Sn(4/ 2 ).l 
shows nearly no increase in M s contradicting the exper- 
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TABLE II. Ten lowest energy configurations of Zn-Sn- 
substituted SrM SrFej^-^Zno.sSno.s^Oig for x = 1. Total 
energies (-Etot) are relative to that of the ground state config- 
uration Zn(4/i)Sn(/2).l. Saturation magnetic moment (M a ) 
and the change in M 3 with respect to pure SrM in Gorter's 
configuration (AM a ) are also given. All values are for a dou- 
ble formula unit cell of SrFei2_a;(Zno. sSno.s^Oig containing 
64 atoms. All energies are given in eV while magnetic mo- 
ments are given in /j,b- 



index 


config 


Etot 


M s 


AM S 


1 


Zn(4/i)Sn(4/ 2 ).l 


0.000 


50.002 


10.001 


2 


Zn(4/i)Sn(4/ 2 ).4 


0.158 


50.001 


10.000 


3 


Zn(12fc)Sn(4/ 2 ).3 


0.177 


40.007 


0.006 


4 


Zn(4/i)Sn(4/ 2 ).3 


0.187 


50.000 


9.999 


5 


Zn(4/i)Sn(12fc).3 


0.190 


40.000 


-0.001 


6 


Zn(4/i)Sn(4/ 2 ).2 


0.257 


50.001 


10.000 


7 


Zn(4/i)Sn(2o).l 


0.265 


39.998 


-0.002 


8 


Zn(2b)Sn(4/ 2 ).l 


0.277 


40.003 


0.002 


9 


Zn(2a)Sn(4/ 2 ).l 


0.277 


40.001 


-0.000 


10 


Zn(12fc)Zn(4/ 2 ).2 


0.384 


40.000 


0.000 



imental result showing a rapid increase of M s with con- 
centration x (~ 7% increase at x = 1). We note that 
configurations with Zn and Sn substituted in 4/1 and 
4/2 sites show considerable increase in M s (~ 10 /xs) 
while those with 2a, 2b, and 12k sites show nearly zero 
increase in AL. 



In Table III we compared the contribution of different 
sublattices to the total magnetic moment in these two 
proposed models for the structure of Zn-Sn-substituted 
SrM. To see the effect of Zn and Sn atoms in differ- 
ent sublattices, we split the entries of sublattices con- 
taining these atoms (4/i, 4/ 2 , and 2b). As expected, 
Zn 2+ and Sn 4+ ions carry negligible magnetic moments 
regardless of their substitution sites. Consequently, 
when they replace Fe 3+ ions in the minority spin sites 
(4/i and 4/2), they eliminate negative contribution and 
make positive contribution. On the other hand, when 
Zn 2+ and Sn 4+ ions replace Fe 3+ ions in the major- 
ity spin site (26), they eliminate positive contribution 
and make negative contribution. Therefore, configura- 
tion Zn(4/i)Sn(4/ 2 ).l results in an increase of 10 \xb for 
the saturation magnetic moment M s , whereas configura- 
tion Zn(26)Sn(4/2)-l, where Zn and Sn are substituted 
to 2b and 4/ 2 sites, shows no increase in M s as the effect 
of Zn in 2b site cancels out that of Sn in 4/2 site. 

In addition to increase of M s , Ghasemi et al. also re- 
ported a rapid decrease in anisotropy of SrM due to Zn 
and Sn substitution. 7 The decrease in anisotropy was at- 
tributed to the substitution of Zn 2+ ions to 2b site. We 
calculated magnetic anisotropy energy (MAE) for the 
two competing models. A self-consistent calculation on 
the completely relaxed crystal structure was performed 



and the charge densities were extracted. Then non-self- 
consistent calculations were performed to structures with 
the spin quantization axis oriented in the easy and hard 
axes. Fig. [5] shows the optimized ground state configu- 
ration Zn(4/i)Sn(4/2)-l with magnetic moments in the 
easy axis and the hard axis. MAE is the difference of 
total energies from these calculations: 



E- 



MAE 



E, 



(100) 



E, 



(001) 



where -E(ioo) is the total energy with spin quantization 
axis in the hard direction (y-axis) and -E(ooi) the total 
energy with spin quantization axis in the easy direction 
(x-axis). Using the MAE, magnetic anisotropy constant 
K can be computed: 



K 



EuAE 

V sin 2 6 



where V is the equilibrium volume and 9 is the angle be- 
tween the two spin quantization axis orientations (90deg 
in this case). Results of these calculations are presented 
in Table |IV| which shows that both the configurations 
show rapid decrease in anisotropy consistent with the ex- 
perimental results. 



IV. CONCLUSION 

Using first-principles total energy calculations 
based on density functional theory, we obtained the 
ground state structure for Zn-Sn-substituted SrM, 
SrFei2-2:(Zno.5Sno. 5)^019 with x = 1. We showed that 
Zn and Sn atoms preferentially occupy 4/i and 4/2 
sites. The stable structure derived from our calculations 
shows a rapid increase in saturation magnetic moment 
and a significant decrease in magnetic anisotropy with 
respect to pure SrM that is consistent with experimental 
observations. We showed that the previously proposed 
model with Zn atom in 2b is an excited configuration 
with significantly higher total energy and consequently 
unlikely to be formed in thermodynamic limit. We also 
showed that the decrease in magnetic anisotropy is not 
associated specifically with 2b site. 
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TABLE III. Contribution of sublattices to the total magnetic moment of Zn-Sn-substituted SrM structures Zn(4/i)Sn(4/2).l 
and Zn(26)Sn(4/2).l compared with pure SrM. All moments are given in fiB- 



site 


SrM 


Zn(4/i)Sn(4/ 2 ).l 


Zn(2fe)Sn(4/ 2 ).l 


atoms 


771 


atoms 




Am 


atoms 


HI 


Am 


2d 




yj .UVJU 


2*Sr 


n no 11 ! 

U.UUt) 




2*Sr 


vj.vjyjo 




2a 


2*Fe 


8.312 


2*Fe 


8.468 


0.156 


2*Fe 


8.313 


0.001 


2b 


l*Fe 


'rk.UO 1 


l*Fc 


4.148 


n dqi 


l*Fe 




u.uuo 


l*Fe 


4 0^7 


l*Fe 




n nm 


l*Zn 


01 ^ 

-vj.vj 10 


4 D70 


4/i 


3*Fe 


-12.114 


3*Fe 


-12.071 


0.043 


3*Fe 


-12.099 


0.015 


l*Fe 


-4.038 


l*Zn 


0.059 


4.097 


l*Fe 


-4.037 


0.001 


4/2 


3*Fe 


-12.285 


3*Fe 


-12.191 


0.094 


3*Fe 


-12.348 


-0.063 


l*Fe 


-4.095 


l*Sn 


0.049 


4.144 


l*Sn 


0.034 


4.129 


\2k 


12*Fe 


50.040 


12*Fe 


50.219 


0.179 


12*Fe 


50.160 


0.120 


4e 


4*0 


1.412 


4*0 


1.251 


-0.161 


4*0 


1.411 


-0.001 


4/ 


4*0 


0.356 


4*0 


0.346 


-0.010 


4*0 


0.333 


-0.023 


6/1 


6*0 


0.164 


6*0 


0.422 


0.258 


6*0 


-0.087 


-0.251 


Yik 


12*0 


1.004 


12*0 


1.034 


0.030 


12*0 


1.016 


0.012 


12k 


12*0 


2.156 


12*0 


3.106 


0.950 


12*0 


2.328 


0.172 


Total 


40.001 




49.957 


9.956 




30.964 


-0.037 



TABLE IV. Magnetic anisotropy energy (MAE) and mag- 
netic anisotropy constant K for two of the x = 1 configura- 
tions of Zn-Sn-substituted SrM and pure SrM. 

Config MAE (mcV) K (kj-m~ 3 ) 

SrM 0.84 190 

Zn(4/i)Sn(4/ 2 ).l 0.45 100 

Zn(26)Sn(/ 2 ).l 0.56 120 
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